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ABSTRACT 

The  weak  form  of  the  free  houndary  problem  for  an  axisymmetric  partially 

penetrating  well  may  be  formulated  as  follows:  find  ^(r)  e  C®([r.,r^])  and 

u  e  C° (ft)  n  V1 (ft)  such  that 

fa  rVu'Vvdrdz  =  0  for  all  v  e  ^ 

and  u  satisfies  appropriate  boundary  conditions .  Here,  u  is  related  to 
the  hydraulic  head,  <fi( r)  is  the  unknown  water-air  interface,  ft  is  the 
region  of  saturated  flow 

ft  =  { (r,z) |0  <  r  <  rQ,  0  <  z  <  h)  U  {(r,z)|r0  <  r  <  r^ ,  0  <  z  <  ^(r)}  , 

K1  is  a  convex  set  in  the  weighted  Sobolev  space  V1 (ft) . 

We  reduce  the  problem  to  three  families  of  variational  inequalities  by 
using  a  type  of  "Baiocchi  transform",  study  equivalence  of  the  three  families 
and  regularity  of  the  solutions  of  the  variational  inequalities.  Finally,  we 
prove  the  existence  of  the  solution  for  the  well  problem. 


AMS  (MOS)  Subject  Classifications:  35J20;  35J65;  35J70;  35R05;  35R35;  76S05. 
Key  Wbrds:  Free  Boundary  problem;  axisymmetric  well;  weighted  Sobolev  spaces; 

families  of  variational  inequalities;  existence. 
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SIGNIFICANCE  AND  EXPLANATION 


When  an  axisymmetric  well  partially  penetrates  a  water  aquifer,  the  water 
flows  through  the  ground  towards  the  well.  By  pumping  water  from  the  well, 
steady  flow  is  obtained.  The  flow  is  governed  by  a  linear  second  order 
elliptic  differential  equation  which  degenerates  at  the  axis  of  symmetry.  We 
reformulate  the  problem  as  families  of  variational  inequalities,  and  study  the 
regularity  of  the  solutions  of  these  variational  inequalities.  Finally  we 
prove  the  existence  of  the  solution  for  the  well  problem.  The  variational 
inequality  formulation  suggests  a  new  numerical  method  for  the  partially 
penetrating  well  problem. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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THE  SOLUTION  OF  THE  FREE  BOUNDARY  PROBLEM 
FOR  AN  AXISYMMETRIC  PARTIALLY  PENETRATING  WELL 


*  •* 

C.  W.  Cryer  and  S.  Z.  Zhou 

Introduction 

The  free  boundary  problem  for  a  fully  penetrating  well  in  a  layer  of  eoil 
of  permeability  K(x,y)  =  exp[f(x)  +  g(y)]  has  been  solved  by  Cryer  and 
Fetter  [1979]  using  variational  inequalities.  In  this  paper  we  consider  the 
problem  for  a  partially  penetrating  well.  A  type  of  "Baiocchi  Transform" 
(Baiocchi  [1974] )  is  used  to  derive  a  corresponding  family  of  variational 
inequalities.  Existence  of  the  solution  is  proved.  To  this  end  we  use  the 
theory  of  weighted  Sobolev  spaces  and  some  results  in  Chang  and  Jiang  [1978]. 

1 .  Weighted  Sobolev  Spaces 

Our  problem  is  governed  by  a  degenerate  elliptic  equation.  Degenerate 

elliptic  equations  can  often  be  associated  with  a  weighted  Sobolev  space  (e.g. 

Murty  and  Stampacchia  [1968],  Trudinger  [1973]).  Various  kinds  of  Sobolev 

spaces  have  been  studied  (e.g.  Jakovlev  [1966),  Cryer  [1980],  Chang  and  Jiang 

[1978],  Leventhal  [1975]  and  Zhou  [1980]).  We  recall  some  results. 

Let  A  be  a  bounded  domain  in  the  (r,z) -plane  with  a  locally  Lipschitz 

00 

boundary  T,  and  with  r  >  0?  CQ(A)  -  the  space  of  functions  infinitely 

00 

differentiable  and  with  support  compact  in  Aj  CgfAjI^)  -  the  space  of 
functions  infinitely  differentiable  in  A  and  vanishing  in  some  neighborhood 
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of  r  #  where  ^  C  r.  L^fAjr)  -  the  space  of  measurable  functions 
satisfying 

Ivl  *  /.  r|v|P  drd z  <  •  . 

LP(A,r)  A 

We  define  weighted  Sobolev  spaces  as  follows: 


(1.1) 


V°(A) 


L2(A;r> 


V^A)  *  {v|3°v  e  L2(A;r),  |a|  <  l) 
v2(a)  «  (v|  ^  32v  e  L2(Ajr),  |a|  <  2) 


(1.2) 


with  norms,  respectively. 


«V  _ 

Ivl 

V°(A) 

LP(A,r) 

Ivl 

1  0 

v\a) 

|a| <1  V°(A) 

Ivl  - 

l  l3°v*  -  + 

V^(A) 

|a|<2  v  { A) 

•  ),  vj( A/J^) 

respectively  the 

1-1  IS, 
r  81  V°(M 


(1.3) 


J 


v’tA). 
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Lemma  1.1.  V  (A),  V  (A)  and  V  (A)  are  Banach  spaces. 

Lemma  1.2.  (Green's  Formula).  If  u,  v  e  V^(A),  then 

t  3v  r  3 ( ru)  r 

JA  ru  drdz  *  v  — —  drdz  +  /j.  ruv  cos(n,r)ds 

r  3 ( ru )  ( 

*  “*a  V  ~3z —  +  Jp  ruv  cos(n,z)ds 

where  n  is  the  outer  normal  of  T, 

Lemma  1.3.  If  A_  is  a  closed  subdomain  of  A  and  3a„  ^  (r  *  0}  *  0, 

—  — 1  -  r— t  fc 


then 


V^Ag)  -  H1(Ae) 
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Now  let  A  be  the  three  dimensional  domain  formed  by  rotating  A  about 


z-axis,  and  let  be  the  surface  formed  by  rotating  about  the  z-axis. 


Lemma  1.4.  If  v(r,z)  e  VK(A),  k  -  0,1,2  and 


/  2  2~ 

f(x,y,z)  =  v(/  x  +y  ,  z) 


(1.4) 


k  *  k.  *  ♦ 

then  f  e  H  (A  ),  where  H  (A  )  is  the  usual  Sobolev  space,  and  A  is  the 

— * 

interior  of  A  . 

Lemma  1.5.  If  v  e  V^(A)  and  (r  =  o)  =  0,  then 

lf«  .  #  -  2TTlvl 

H  (A  )  V  '  (A) 

/s  f2dS  *  2 w  /p  rv2ds  . 

By  using  Lemma  1.5  and  results  in  Sobolev  [§10,  1950]  we  obtain: 

Lemma  1.6.  If  v  e  )  and 

mestrN^r.  n  {r  =  o})]  >  0 


then 


>V'2,  <  C  /a((£)2  .  (£)2,r  dr*. 

V  (A) 


where  C  does  not  depend  on  v. 


2.  Desoritption  of  the  Problem 

The  problem  to  be  considered  is  shown  in  Figure  2.1. 

A  cylindrical  well  of  radius  rQ  partially  penetrates  a  layer  of  soil  of 
depth  H  and  radius  r^  Take  the  axis  of  symmetry  as  the  z-axis.  The 
bottom  of  the  soil  layer  is  impermeable.  The  distance  of  the  well  bottom  from 
the  bottom  of  the  soil  layer  is  h.  We  assume  that  the  soil  layer  is 
homogeneous  and  isotropic?  that  the  water  is  imcompressible?  that  the  flow  is 
irrotational  and  steady  (in  particular  the  height  of  water  on  the  outer 
boundary  of  the  soil  and  in  the  well  is  respectively  H  and  hw);  that  the 


permeability  k(r,z)  =  1 


and  z  «  <fi(r)  ia  the  boundary  between  the  wet  region  ft  *  ft^  U  ft^  and  dry 
region  ft3»  It  is  called  the  free  boundary  as  it  is  the  unknown  part  of  3ft. 

Denote  by  p(r,z)  and  u(r,z)  respectively  the  pressure  at  point  (r,z) 
of  D  (the  atmospheric  pressure  being  measured  by  zero)  and  the  hydraulic 
head,  then  we  have 

u(r,z)  *  p(r,z)  +  z  in  ft  .  (2.2) 

It  follows  from  Darcy's  law  and  the  equation  of  continuity  that  (see  Hantush 
[1964],  or  Cryer  (1976,  p.  86]) 

Lu  ,  ♦  i  ♦  fa .  „  ln  a  .  ,2.3) 

3r2  r  3r  3z2 

We  introduce  the  notation 

r  *  {(r,z) | 0  <  r  <  r^,  z  -  0} 

I*  »  {(r,z)|r  *  r1#  0  <  z  <  H> 

r3  -  {(r,z)|rQ  <  r  <  rt,  z  «  H} 

T  ■  {(r,z)|r  -  r  ,  h  <  z  <  H> 

4  0  w 

r  =  {(r,z)|r  «  r  ,  h  <  z  <  h  } 

3  0  W 

T  =  {(r,z)|0  <  r  <  r  ,  z  =  h) 

6  0 

r?  =  t(r,z)|r  »  0,  0  <  z  <  h} 

»  {(r,z)|r0  < .  r  <  r1 ,  z  =  v’(r)}. 

Then  u(r,z)  satisfies  the  following  boundary  conditions: 


u  * 

H 

on 

r2 

(constant  hydraulic  head)  > 

u  * 

z 

on 

r  u 
0 

( r  o  3ft) 

4 

(interface  with  air) 

^  (2.4) 

u  * 

h 

w 

on 

r  u 

5 

r 

6 

(interface  with  water  at  rest) 

n 


V  i 


3u 


t—  ■  0  on  r  U  T  (streamline) 
on  0  1 


3u 


(2.5) 


0  on  r. 


(symmetry) 


Now  we  can  state  our  problem  in  weak  form. 

Problem  (PPW) 

Given  the  domain  as  in  (2.1),  and  a  real  number  hw  such  that 
h  <  hw  <  H,  find  functions  <i(r)  and  u(r,z)  such  that  u  satisfies  (2.4) 
and 


<fi  e  C  (Ir0,r1l),  (r^ )  *  H,  so  ( rQ )  >  hw 

<P  is  strictly  increasing  , 
u  e  v^Q)  n  c°(T5)  , 

/jj  rVu*Vv  drdz  *  0  for  all  v  e 


(2.6) 

(2.7) 

(2.8) 
(2.9) 


where 


fl  *  fli  u  {(r,z)  e  D|r  >  rQ,  0  <  z  <  ^(r)} 

K  -  (v  e  vVfl^v  =  0  on  r  u  ( T  n  3(2)  u  V  r  }  . 

i  2  4  5  6 

Remark  2.1.  This  problem  can  be  regarded  as  a  plane  problem  with  permeability 
K  =•  exp ( in  r) ,  but  it  is  not  covered  by  the  work  of  Benci  [1974]  because 
In  x  $  * 2+M( [Q,r^] ).  Also,  it  can  not  be  included  in  Alt  [1979]  as  a  two- 

dimensional  problem.  Rama  and  Das  [1976]  solved  this  problem  by  numerical 
methods. 

Remark  2.2.  Chang  and  Jiang  [1978]  have  solved  a  similar  problem  by  using  so- 
called  "Sequence  of  set-valued  mappings"  instead  of  the  method  of  variational 
inequalities.  Further  results  about  obstacle  problems  have  been  obtained 
recently  by  Chang  [1980].  We  use  some  results  on  linear  equations  given  by 
Chang  and  Jiang  [1978].  But  we  solve  our  problem  by  using  the  method  of 
variational  inequalities  because  the  corresponding  numerical  method  is  more 
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convenient,  and  because  in  our  case  the  boundary  conditions  and  right  term  of 
the  nonlinear  equation  for  the  Baiocchi  function  w  are  different  from  those 
in  Chang  and  Jiang  [1978]. 


3.  The  Baiocchi  Function  w  and  its  Properties 

Assuming  a  priori  the  existence  of  the  solution  u  of  (PPW),  set  (see 
Baiocchi  [1974],  [1976],  [1978]] 


:,z)  =  ^ 


u(r,z)  in  A 
z  in  D\fi 


(3.1) 


w(r,z)  =  /*  [u(r,t]  -  t] dt 


(3.2) 


Remark  3.1.  We  can  not  use  the  transform 


fv’l(r)  - 

w(r,z)  =  J  (u(r,t)  -  t] dt  , 


where  ^(r)  =  H  for  rQ  <  r  <  r1  and  ^(r)  =  h  for  0  <  r  <  rQ,  since 
3w  t  1  — 

^  f  C  (D).  This  is  obvious  physically.  (Cf.  Lemma  3.6). 

Now  we  derive  some  properties  of  w  and  u. 


Lemma  3.1. 


Lu  =  0  in  0 


(3.3) 


u  is  analytic  in  8\{r  ,  a,  b,  C,  F,  G,  P,  q) 


(3.4) 


g  -  0  on  r,  u  r7 


(3.5) 


x—  =  0  in  the  weak  sense  on  I  . 
3n  0 


(3.6) 


Proof .  On  writing  (2.9)  for  any  v  e  Cg(fi),  we  obtain  (3.3)  in  the  sense  of 
distributions.  From  classical  results  on  the  regularity  of  the  variational 
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solutions  of  elliptic  equations  in  the  interior  and  on  the  smooth  parts  of  the 


boundary  (see  for  instance  Lions  and  Magenes  [1972,  §9,  Ch.  2])  it  follows 
that 

u  is  analytic  in  ^  *  A,  B,  C,  F,  G,  P,  Q }  . 

— * 

Denote  by  ft  the  three-dimensional  axisymmetric  domain  with  cross-section 

7%  *  /  2  2 

«•  Then  u  (x,y,z)  =  u(»  x  +  y  ,  z)  is  a  solution  of  the  equation 

★  *  _ *  £ 

Au  *  0  in  ft  where  8  is  the  interior  of  ft  .  Hence  u  and  u  are 

* 

analytic  on  because  1*^  is  in  ft  ,  and  (3.4)  is  valid.  (2.9)  implies 

(3.5)  and  (3.6)  in  the  weak  sense.  (3.3),  (3.5)  are  satisfied  also  in  the 
classical  sense  thanks  to  (3.4). 

Q*  E • D • 

Remark  3«2«  It  follows  from  the  three-dimensional  argument  above  and  maximum 

principle  that  if  Lv  >  0  in  ft  with  v  e  C°(ft)  then 

v| -  <  max  v 

7  ft 

3  v 

and  if  in  addition  <  0  at  the  point  Q  then 

v|  <  max  v 

Q  a 


Similar  results  are 

valid  for  min 

V 

if  Lv  >  1 

Lemma  3.2. 

u(r,z) 

>  2 

in  ft  . 

Proof.  Set  v  =  u  - 

z .  Then  we 

have 

r 

Lv  *  0 

in 

ft 

v  =  0 

on 

ro  Oft 

J 

v  =  H  -  z 

on 

r 

2 

\ 

v  =  h  -  z 
w 

on 

r  u  r 

5  6 

3 v 

1  ^  =  +1 

on 

ri 

3v 

^  3n  =  ° 

on 

r7  * 

(3.7) 
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Since  ti  is  elliptic  and  v  e  C^(ft),  v  attains  its  minimum,  m  say,  in  ft 

*  *  .  *  / 

at  a  point  p  e  oft.  But  p  f.  (by  Hopf  principle)  and  p  $  (by 

Remark  3.2).  So  p  e  3ft\(T^  ) .  Hence  m  is  zero.  It  follows  from  the 

strong  maximum  principle  that  v  >  0  (i.e.  u  >  z)  in  ft. 

QiG»D* 

Lemma  3.3.  u  8  v’(D)  n  C°(D)  (3.8) 

-  3*ft 

m  »  -  t —  in  the  sense  of  distributions  (3.9) 

®z 


where  is  the  characteristic  function  of  ft  in  D. 

Proof.  By  (2.4),  (2.8)  and  (3.1)  it  is  easy  to  see  that  u  e  C°(D).  For  any 
00 

4*  e  Cg(D)  we  have 

4>  "  I!  drdz  =  ^  +  + 

=  C  dr  fo  u  H  dz  +  dr  Cr)  u  H  dz  +  dr  ^(r)  a  ^ dz 

9\i  3vi 

=  -/j j  41  drdz  -  /jj  4>  jp  drdz  -  Jq  V  drdz  (Integration  by  parts) 

1  Z  2  Z  3 

=  /D  t  v  drdz 

where 


in  ft 


Hence 

in  ft 
in  D\ft 

Similarly  we  have 
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^  0  in  D\n 


Clearly,  u  e  V  (D). 

00 

Now  we  prove  (3.9).  If  ^  e  CQ(D)  then  <|>  e  .  Hence  we  have  for  the 

— 00 
distribution  Lu  and  every  e  Cn(D) 

<Lu,r<|»>  ■  -/  rVu*Vi|>  drdz 


-  fa  rVu'V^  drdz  -  ||  r  drdz 

/  3ili  a  3 

-  4  ^ (1  ■  Vr  ^  (1  ■  V'  r*> 


=  -  <-5 — ,  ri l»  . 


It  is  just  (3.9). 


Q«E«D» 


Proposition  3# 4%  Let  w  he  defined  by  (3#2)*  Then 

Lw  =  in  the  sense  of  distributions. 

Proof .  Since  u  e  C°(D),  we  have 

3w  * 

*  u(r,z)  -  z  . 

So  in  the  distribution  sense  we  obtain  (by  (3.9)) 


(3.10) 


(3.11) 


3  ^3wi  “ 

«-  (Lw)  -  =  Lu  -  Lz  =  - 


3z  * 


Hence  (Schwartz  [1973,  §5,  Ch.  2]) 


Lw  +  *Q  =  T(r)  ft  1(z)  . 


Since  u  is  analytic  in  n  u  r  and  *  0  on  we  have  in  Os 


3^w  1  3w  3^w 

Lw  *  - „  +  -  ■5—  +  — - 

Sr2  r  3t  3z2 


/;  (La  ♦  i  |a)dt  *  Lit, 

0  r  3r  3z 

rz  _  ..  3u(r,0) 

*  I.  Lu  dt  +  — f  |(|  ■  -  1  ■ 


Accordingly,  T(  r )  =  0  and  (3.10)  is  valid.  . 


Q.E.D. 


Proposition  3.5.  w  is  a  solution  of  the  equation  with  discontinuous 
nonlinearitiess 


Lw  * 


for 

for 


3w 

3z 

3w 

3z 


-  0 
>  0 


Proof .  It  is  enough  to  prove  that 

If  >  0  '  lf  <r'z)  e  °  ' 

A 

2^  ■  0  ,  if  (r,z)  e  d\8  . 

But  this  is  obvious  thanks  to  Lemma  3.2  and  (3.11). 


(3.12) 


Lemma  3.6. 


3w 

37 


tt  3u 
0  3r 


dt 


Q.E .D. 


(3.13) 


r  *  constant  *  0  for  r  e  Ir^r^  .  (3.14) 

00 

Proof.  For  any  e  CQ(D)  we  have 

ft ^ w  If drdz  =■  /J dz  f0° [*(r,t)  ■ t]  If dr 

=  /;  dz  /*  dt  /0°  (u(r,t)  -  t}  If  dr 

*  J|J  dz  JJ  [u(r,t)  -  tl ♦!  dt  -  /a  *</J  |-  dt ) drdz 

0  1 

(integration  by  parts). 


Similarly,  we  have 


W  U  drdZ  "  dZ  ^  "  t]Mr=r0At  -  W,  K  S  dt)drd2 


3u 


Since  ^  *  0  on  {(r#z)|r»r^#  h  <  z  <  h),  we  obtain  immediately 


/p  W  ^r<*z  m  ~/D  ♦</?  — <3t)drdz  for  any  +  6  Cg(D) 


^w( r  H) 

(3.13)  has  been  proved.  Now  set  f(r)  =*  r  — jj— * — Then  for 


any 


'Hr)  e  C  Or  ,r  [>  we  have 


V  If  <r>f(r>dr  =  /j.1  If  r  dr 

0  0 

-  i! 

-  £  *•  /;  g  *  «■ 

0 

=  /n  r^u’V^  drdz 


where 


^(r.z) 


{ 


<|»(r)  for  r  e  ]rQ,r1  [ 


for  r  e  ]0,  r0) 


Clearly  ^  e  K^,  Hence  it  follows  from  (2.9)  that 


J  1  f  dr  *  0  V^e  C  Or  ,r  I) 
r^  3r  001 


so  f(r)  is  a  constant  for  r  e  Jr^r^!,  which  we  denote  by  q. 


Q.E 


Remark  3.3.  Physically,  2nq  is  the  discharge  of  the  well. 
Proposition  3.7.  Let 
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0 


Hz  -  T 


H  »  r 

—  +  q  *n  — 
2  r . 


on  r  u  r 

3  4 


(3.15) 


H2  t  .  'O  (VZ)  - 

r  t  ,  in  - - 5 -  on  r5 


h  -  h  on  T 
w  6 


on  r 


(3.16) 


w(r,z)  =  g  on  T 

q  D 


3n  % 


where  ^  ^  1^,  ^  =  Tg  U  r?. 

Proof.  Thanks  to  (3.2)  and  (2.4)  we  have  clearly 


w  =  0  on  ri#  w  -  H*  -  —  on  I* 


Hence  w(r.,,H)  =  H2/2.  Solving  the  ordinary  differential  equation 


3w(r ,H)  .  .  . 

r  — 1^7* —  =  q  we  obtain 


H  .  r  r 

w  =  —  +  q  *n  —  on  I 
2  ^  r  3 


— • 

On  r  we  have  -s—  =  u  -  z  =  0.  Hence 
4  °z 


2  r 

w(rQ,z)  =  w(rQ,H)  =  j-  +  q  *n  ~  . 


On  T  we  have  t—  *  h  -  z  and 
5  oz  w 


,,  3w  H2  r0  (V2) 

w  *  w(rQ,hw)  +  Jh  ^  dz  -  j-  +  q  *n  - - - 

w  1 


-13 


(3*16)  is  obvious. 


Proposition  3.8. 


w(r,z)  =  g  (r,H)  in  D\ft 

q 


w(r,z)  <  g  (r,H)  in  ft\ft,  . 

q  1 


(3.17) 


(3.18) 


Proof .  (3.17)  follows  from  (3.1),  (3.2)  and  Proposition  3.7.  (3.18)  follows 


from  (3.17),  Lemma  3.2  and  the  fact  that  u  B  C  (ft). 


QsEeD* 


Remark  3»4,  We  obtain  another  form  of  the  nonlinear  equation  for  w: 


0  in  {w  *  g  (r,H)} 

q 


Proposition  3.9. 


-1  in  ft,  U  (w  <  g  (r,H)} 

'  q 

w  e  v2(d)  . 


Proof.  By  (3.2)  and  (2.8)  we  have 


3w  d2w  32w  0.  . 

K'  £5'  e  V  (D) 


(3.20) 


(3.21) 


Differentiating  (3.13)  we  obtain 


Now  we  prove  that 


In  fact,  we  have 


3  w  3u  „0/P.v 

Ttt  -  7?  e  V  (D)  * 


3w  0 

^ev°<D)  . 


(3.22) 


(3.23) 


/  r(/ *  dt)2  drdz  <  /n  rH [/*  (t^)2  dtjdrdz  (Schwartz  inequality) 


'D  L"l'0  '•3r 


"i/0° dr  4  *  /?  *  *  4’ dr  /; dz  4 r^2  at> 


<  HI-  4°  ar  4  r(f»)’«  *  «  4’  d-  4  r(g)2  *, 

o 


<  H2  /D  r(|^)2  drdz  <  •  . 


At  last  we  prove  that 


1  3w  3_w  0,  . 

tj—  t  ~  6  V  (D)  • 

r  or 


(3.24) 


By  (3.10)  we  have,  as  distributions, 


32w  »  32w  1  3w 

3r2“'n"3z2"?ar  ' 


Hence  it  is  sufficient  to  prove  that 


1  3w  _  „0 
-  ^  e  V  (D) 


(3.25) 


Let  R*  =  C(r,z)|0  <  r  <  0  <  z  <  hi,  Q*  «  <~,h),  A*  -  (~,0),  f(z) 


u(—,z).  Then  u|_^  e  C  (R  \(a,Q}>  n  c  (R  )  n  v  (R  )  (by  Lemma  3.1),  and 
2  R 

u|  *  is  the  solution  of  the  boundary  value  problem 
R 


Lu  =  0 


* 

in  R* 


u|  ,  =  h  ,  -s—  |  .  =  0 

z=h  w  3z  z=0 


u|  =  f(z),  u  is  bounded  near  r  *  0  . 

ro 

r~r 


By  using  the  method  of  separating  variables  we  obtain  that 


V«>  - .  I  ilr  tj4ht  (f)*  • 


(3.26) 


3n  1 

It  is  easy  to  show  that  ^  *  0(r>  as  r  +  0*  Hence  -  *  0(1)  as 

r  *  0,  and  -  .  e  V°(R  ).  How  (3.25)  is  clear, 

r  3r  5* 


Q*  E • D« 
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Proposition  3.10. 


w  e  C^D)  . 


(3.27) 


In  order  to  prove  this  proposition  we  cite  a  theorem  in  Chang  and  Jiang 

11978],  the  proof  of  which  is  given  in  the  Appendix. 

Theorem  3. A.  If  v  e  V2(D)  and  Lv  =  f  e  LP(Djr)  (p  >  6), 

v|r  =  ^Ir  “  then  v  e  C^(D)  (0  <  \). 

D  3n  N  2 

Proof  of  Proposition  3.10: 

At  first  we  construct  a  function  v  such  that 

q 


v  lr  = 
q  d  q 


v  e  v2(d)  n  c1  (d) ,  lv  e  L^D/r)  . 
q  q  ' 


(3.28) 


To  this  end  we  set  (cf.  Baiocchi  et  al.  [1973],  Chang  and  Jiang  [1978]) 


v  «  vi  +  qv2 


(3.29) 


where 


f0(z) 


2r-r  -r  2 

[f  (z)  -  f0(z)]  (— — l — -)  +  fQ(z) 

10 


0  <  r  < 


ro+ri 


r  +r„ 

0  1  „ 

-  s  r  S  r 

2  1 


(3.30) 


f2(r)  [1  -  (~)2J  0  <  z  <  h 


f2(r) 


h  <  z  <  H 


(3.31) 


f(z)  = 


(h  -z) 

w 


h  <  z  <  H 
w 


0  <  z  <  h 
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f(z) 


h  <  z  <  H 


Vz)  = 


f(z)[1  -  (^)21  o  <  z  <  h 


f^z)  =  HZ  -  ~ 


f2(r)  = 


'  fa  ^  ♦  3(^-)2 
ri  0 


in  — 
r1 


2(H3 

ro 


-  1  0  <  r  <  rr 


ro  <  r  ‘  ri 


(3.32) 


It  is  readily  verified  that 


vi'r 


»*  -! 


H 

2 


2  (h  -z) 
H  w 


k  2 


on 


on 


on  r3ur4 


on 


(3.33) 


9v 

*T'r_  '  %> 


N 


(3.34) 


v2lr 


^  ri 


on  rn  u  r2 


in  z-  on  r3  U  r4  u  r 


(3.35) 


3V2 

■jr-l  r  “  0 

3n  N 


(3.36) 


v  ,  v  e  v2(d)  n  c\d ) 
1  2 


(3.37) 


Lv^ ,  Lv2  e  L  (0;r)  . 


(3.38) 
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Hence  (3.28)  is  valid 


Now  we  set  v  =  w  -  vq,  then  v  e  \r<D)  (by  (3.20)  and  (3.28)),  and 

Lv  =  -  l>v^  e  L  (Djr) 


v,r  =  fe'r  =  0 

D  N 


(3.27)  follows  from  Theorem  3. A  immediately. 


Q.  E  »D  . 


Proposition  3.11 


q  =*  lim  — /”  Iu(r,t)  -  u(r  ,t))dt 
j.  ,  „  r-r  n  0 

r+rQ+0  0 


+  lim  — ~  /  °[u(r,z)  -  u(r,h)]r  dr  . 
„  z-n  u 
z+h-0 


(3.39) 


Proof :  By  the  mean  value  theorem  of  differentiation  we  have  that  for  any 
r  e  ]rQ,r  [  there  exists  £  e  ]rQ,r{  such  that 


Jh  Iu< r ,t)  -  u(rQ,t)]dt  -  /h  [u(r,t)  -  tjdt  -  /htu(rQ,t)  -  tjdt 


w  ( r , H)  -  w(r,h)  -  [w(rQ,H)  -  w(rQ,h)J 


3w(g,H)  3w(C,h) 

1  3r  3?  J  U  "  V  * 


It  follows  from  (3.27)  and  (3.15)  that 


lim  — /”  [u(r ,t) 
r>rQ+0  r"r0  h 


u(r0,t))dt  *  q  -  rQ 


3w(rQ,h) 

37 


(3.40) 


On  the  other  hand,  for  any  z  e  ]0,h[  there  exists  a  n  e  ]z,h[  such  that 


/0°[u(r,z)  -  u(r,h)]r  dr  *  /Q°  r  dr  . 


(3.41) 
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Let  D1  «*  {(r,z)|0  <  r  <  rQ,  0  <  z  <  n).  It  is  easy  to  show  by  Green's 
formula  that 


/  Lu*r  drdz 


r 

'o 


3u(rQ,z) 

3r 


dz 


-/r° 

;0 


c2»l £jSIdr 


So  we  have 


rro  hSjLiHl 

0  3z 


r  dr  =  r 


0 


3u(r  ,z) 
fn  0 

J0  3r 


dz 


3w(rQ,h) 
r0  3r 


(3.42) 


Now  we  obtain  (3.39)  by  (3.40)  -  (3.42)  and  (3.27). 


Qe  E* D • 


Remark  3.5.  Physically,  (3.39)  means  that  the  total  discharge  to  the  well 
consists  of  two  parts:  one  is  across  the  wall  of  the  well,  another  is  across 
the  bottom. 


4.  Variational  Inequalities  (VI)  Satisfied  by  w;  Regularity  of  the  Solution 
of  VI 

Define  functions  for  every  v  e  V^(D)s 


v'  = 


V 

9q(r,H) 


yw  a 


C 

I  V  -  g  (r,H) 

'  g 


in  0,  (v  <  gqH> 
in  (v  >  gqH} 
in  Q,  Cv  <  gqH} 
in  {v  >  gqH> 


where  (v  <  g  )  =  (<r,z)  e  D|r  >  r  ,  v(r,z)  <  g  (r,H)} 
qH  u  4 

(v  >  g  }  =  {(r,z)  e  D|r  >  r  ,  v(r,z)  >  g  (r,H)J  . 
qH  0  q 

Then,  clearly,  we  have 


(4.1) 


(4.2) 
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I 


V  -'-/'-"'Jf 


0 

v  =  v*  +  v",  v"  >  0,  v'  <  g  (r,H)  for  r  e  )r  ,r  [,  v',  v"  e  V  (D)  . 

q  0  1  <4. 3) 

Let 

K  =  (v  6  V1(D)jv  =  g  on  T  }  .  (4.4) 

q  q  D 

We  have 

Theorem  4.1.  If  u  is  a  solution  of  (PPW),  then  w  defined  by  (3.2)  is  a 
solution  of  the  VI: 

w  e  K 

q 

r 

/  rVw*V(v-w)drdz  -  (h  -h)  /.  r(v-vi)|  .dr  -  /  (v'-w*  )r  dr  dz  >  0(4.5) 

O  w  0  z-h  D 

for  veK 

q 

Proof :  By  (3.20)  we  have  w  e  V2(D).  Apply  to  w  and  to  any  v  e  Kg 
the  following  Green's  formula: 

/  r^w*V(v-w)drdz 
D 

=  -  /D  rLw*(v-w)drdz  +  jp  r(v-w)  ds  +  jp  r(v-w)  ds 

r 

=  f  r$0* (v-w)drdz  +  (h  -h)  /  °(v-w)|  r  dr 
D  “  w  o  z=n 

r 

«  /a  r(v-w)drdz  +  (h^-h)  JQ  (v-wJl^r  dr 

r  ,r0  .  /-v"  ^  0  and> 

>  r(v'-w')drdz  +  (h^h)  JQ  (v-w)lz=hr  dr  =  p  J  . 

But  -  JD  ^(v'-wMr  drdz  =  -  /D^[v'-g^(r,H)]r  drdz  >  0  (by  (3.17)).  Hence 
(4.5)  is  valid. 

Q*  E  •  D« 
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Remark  4.1.  It  is  easily  seen  (by  (3.17)  and  (3.18))  that  w  is  also  a 
solution  of  the  VI 


v  e  K 


/  r^w*?(v-w)drdz 


(h  -h)  f  (v-w)|  r  dr  -  f  r(v-w)drdz  >  0  (4.6) 

w  0  z*h  D 


for  v  e  K 


where 


C  ■  (v  e  V1  (D)  |v  =*  g  on  r  ,  v  <  g  (r,H)  in  d\0.)  .  (4.7) 

q  q  D  q  1 


q  'q  D  'q 

Remark  4.2.  Noting  (3.7)  we  have  that  w  is  also  a  solution  of  VI 


w  e  K 


r 

/  rVw*V(v-w)drdz  -  (h  -h)/_°(v-w)|  r  dr  -  /  (v-w)r  drdz  >  0  (4.8) 

D  W  0  Z*n  D 


for  vex 


where 


K**  «  (v  e  K*|v  >  0  in  d) 
q  q 


(4.9) 


Remark  4.3.  For  numerical  solutions  (4.8)  is  the  most  convenient  VI.  By  the 


well-known  result  (Lions  [1971,  p.  9]),  Problem  (4.8)  reduces  a  minimization 


** 

problem  on  a  convex  set  as  follows:  find  w  e  such  that 

J(w)  »  min  J(v) 


vex 


** 


(4.8*) 


where 


J(v)  “  f  rl^vl^drdz  -  (h  -h)  v  r  dr  -  2  f  vr  drdz 
D  w  0  z*h  D 


(4.8*)  is  the  basis  of  numerical  solutins  to  (PPW)  by  using  Vi's. 
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For  q  e  R,  (4.4)  is  a  family  of  Vi's 


So  are  (4.6)  and  (4.8).  Now  we 


study  these  families. 

Proposition  4.2.  V  q  e  R,  (4.4)  has  unique  solution  W^. 

Proof:  Set  V  =  (v  e  v\d)|v  =  0  on  r  }  =  V*(D»r  ) 

a(u,v)  =  /  r^u*Vv  drdz 

r 

f(v)  *  /  v’r  drdz  +  <h  -h)  /  rv|  dr  . 

D  w  u  z=n 

Then  V  is  a  Hilbert  space  with  inner  product  (u,v)y 

=  /  r(uv  +  Vu»Vv)dr  dz;  K  is  a  closed,  convex,  non-empty  (e.g.  v  e  K ; 
d  'l  q  g 

see  (3.28))  subset  of  V;  a(u,v)  is  a  bilinear,  continuous  and  coercive  form 
on  V  *  V  (by  Lemma  1.6);  and  it  is  easy  to  show  that  f(v)  is  a  convex, 
continuus  functional  on  V  with  f(v)  ¥  -°®  and  f(v)  ¥  +0D*  By  the  well- 
known  theorem  (Lions  and  Stampacchia  [1967,  theorem  2.2)),  we  obtain  the 
conclusion  of  our  proposition. 

Q.E.D. 

Proposition  4.3.  V  q  <  ,  where 


2  2 
H  -  (h  -hr 

_ w _ 

q0  ~  2  Hn(r1/r()) 


(4.10) 


(4.6)  has  a  unique  solution.  (4.8)  also  has  a  unique  solution. 

The  proof  is  similar  to  that  of  Proposition  4.3.  The  condition  q  <  qQ 

*  •* 

ensures  that  both  K  and  K  are  non-empty. 

q  q 

We  will  prove  later  that  the  problems  (4.4),  (4.6)  and  (4.8)  are 
equivalent  for  q  <  q^  (Theorem  4.13).  Now  we  study  (4.4)  in  detail. 
Proposition  4.4.  V  q  e  R,  the  solution  w^  of  (4.4)  satisfies,  in  the  sense 
of  distributions,  that 


- 1  <  Lw  4  0 

q 


(4.11) 
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I 


& 


i 


V '  -  w*  = 

q 


Hence 


Lw  e  L 
q 

(D;r) 

• 

(D),  tp  >  0.  Let 

v  = 

w  -  tp. 

q 

in 

n  u 

i 

W  -  ♦  -  g  ( r ,  H ) 

q  q 

in 

{9qH 

0 

in 

(w  > 

q 

v'-w* 

q 

>  <P 

• 

v  =  Wg  -  <p  and  w  = 

we 

(4.12) 


•>  s  g  „•< 
q  qH 


q  qH 

3  +  ii>) 

*qH  y 


(4.13) 


r 

0  <  “/n  rVw„ *V^  drdz  -  (h  -h)  /  0  r+|  .dr  -  f  r(v’-w')drdz 
u  q  wo  z=h  D  q 

<  -/_  rVw  *Rip  dr  dz  +  /  r^  drdz  =  <Lw  +  1,  ri|» 

D  q  D  q 


for  any  ♦  e  CQ(D),  \p  >  0  . 


Hence 


Lw  +  1  >  0  . 

q 


(4.14) 


Similarly,  given  ♦  e  C.(D),  ip  *  0,  let  v  =  w  +  ip.  Then  v'  >  w' ,  and 

0  q  q 

(4.5)  becomes  that 


0  <  /  rVw  ’Vipdrdz  -  /  r( v'-w*  )drdz  <  /  rVw  ‘Vipdrdz 
D  q  D  q  D  q 


<Lw^,  r<P  >  <  0  for  any  tp  e  CQ(D),  ip  >  0 


Hence  Lw  4  0. 

q 

This  inequality  and  (4.14)  prove  (4.11),  and  (4. 12 )follows  from  a  well-known 
theorem  (Schwartz  (1973,  th.  v,  p.  29)  and  Radon-Nikodyn  theorem. 

Q.E.D. 


roposition  4.5.  If  wq  is  a  solution  of  (4.4),  then 


3n  *r  "  gN  ' 
N 


Proof:  At  first  we  prove  that,  in  the  sense  of  distributions. 


(4.15) 


v  ■  ^  =  h  -  h  on  T 
az  w  6 


"3?  =  0  on  r7  * 


(4.16) 


(4.17) 


Given  e  C  (T  )  such  that  <J»  >  0,  and  e  >  0,  we  construct  an  element 
0  6 

4»e  e  v\d)  with  4»e  >  o  in  D,  (|>£  =  <|>  on  V&,  i|»e  =  0  on  and 


/  r^  drdz  <  E  . 
O  £ 


(4.18) 


Let  v  =  w  -  *l>  .  Then  we  have  v'  -  w'  >  ♦  similar  to  (4.13).  It  follows 
q  e  q  e 

from  (4.5)  and  generalized  Green's  formula  (Baiocchi  and  Capelo  [1978; 
Appendix  4  of  V. 1 ) 

r 

0  <  -fo  rVw  'Vif^drdz  +  <hw~h)  fQ°  rl^elz=hdr  +  JD  r^e  drdz 


<  <LW  ,n|»e>  -  <-§-^a,r^e>r  +  <\~h‘  r^e>r  +  <1»  1'*e> 


vw 

<x-a  -  (h  -h),  r<|»  >r  <  <LW  +  1,  r  ♦  > 
az  w  E  ‘ ,  q  E 


On  the  other  hand,  writing  (4.4)  with  v  =  wq  +  we  obtain 

r 

0  <  /  rVw  *Vi|»  drdz  -  (h  -h)  /  r4»  I  dr  (since  v'  >  w’) 

'n  q  £  w  0  e  z=n  q 


(4.19) 


-<Lwq'rV  +  <3 jf.  r*e>f 


<h  -h,  ri|»  >- 
w  £  :  , 


i«e. 


<rsf-  <Vh>'  rVr  >  <Lv  rV  • 

6 


(4.20) 


By  (4.19),  (4.20),  (4.11)  and  (4.18)  we  have  (since  i|>  *  ib  on  T  ) 

E  D 


I  3w  1 

"  (hw"h)'  I  <  4  r^s  6x62  <  e  ' 


Since  £  >  0  is  arbitrary,  we  have 


<jf  -  <Vh)*  rVr  =  0  v  ♦  e  co(V'  *  *  0 

6 


This  proves  (4.16).  Now  we  prove  (4.17). 

Introduce  the  notation: 

F  =  Lw 

q  q/ - 

*  /  2  2 

v  *  v(»  x  +y  ,z)  for  any  function  v(r,z) 

*  /  2  2 
D  =  {(x,y,z) | ( r,z)  e  D,  r  =  /  x  +y  }  {(x,y,z) |x=y=0,  0  <  z  <  h} 

=  {(x,y,z) | (r,z)  e  1*^,  r  =  /  x2+y2} 

rD  =  t(x,y,z) | (r ,z)  e  r  ,  r  =  J x2+y2}  . 

* 

Then  D  is  a  three-dimensional  axisymmetric  domain  whose  boundary  is 

*  *  * 

r  U  r  and  w  is  the  solution  of  the  problem: 

D  6  q  r 


Aw  *  F 

q  q 

in 

D 

*  * 

w  =  g 

q  q 

on 

* 

r 

D 

* 

3w 

r-3  «  h  -  h 
3n  w 

on 

r6 
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By  (4.11),  F  6  L  (D  ).  Hence  w  |  e  H^'^(D.).  p  <  “,  where  D. 

q  *  /  2 — 2  q  Dl  1 

{(x,y,z)  |  (x,y,z)  e  D  ,  /  x2+y2  <  rQ  -  $}. 


w*  e  c^d, ) 
q  i 


By  embedding  theorem  we  have  that 

(4.21) 


Now  it  is  easily  seen  that 


*  * 

3w  3w 

3^  *  9^  3  0  at  x 


0  in  D . 


Hence 


3w 

1”^  *  0  at  r  =  0  in  D 


(4.17)  has  been  proved.  Moreover,  (4.21)  means  that  (4.15)  is  valid  in 
ordinary  sense. 

Q.E.D. 

Now  we  need  the  following  results  (see  Chang  and  Jiang  [1978]). 

Lemma  4. A.  Let  f  e  LP(Djr),  p  *  2.  Then  the  problem 


Lv  =  f  in  D 

v’r  "  3n  r  =  0 

D  N 

has  unique  weak  solution  v  in  V^D),  and  v  e  C°(D).  Moreover,  the  linear 
operator  K:  f  v,  mapping  Lp(D;r)  (p  >  2)  into  C°(D),  is  compact. 
Theorem  4.B.  Let 


a  a 

U(D)  *  (v  e  V  (D) | v  =0  on  r  ■—  =  o  on  r  }  , 

u  on  n 


Denote  by  R(L)  the  range  of  the  operator  L  as  a  map  from  U(D)  into 
V°(D).  Denote  by  R(L) ^  the  orthocomplement  of  R(L)  in  V°(D).  Then 

dim  R(L)1  =*  1 
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i.e.  there  exists  vQ  e  R(L)  such  that 

r(l)1  =  {v  e  v°(d)|v  *  wvQ ,  v  e  rI  . 

Remark  4.4.  It  is  easy  to  show  that  K  is  also  a  compact  operator  mapping 
into  V^D). 

Now  we  verify  the  continuity  of  the  solution  of  (4.4). 

Proposition  4.6.  If  wq  is  a  solution  of  (4.4),  then 

w  e  C°(D)  .  (4.22) 

q 

Proof ;  Let  f  *  Lwq  -Lvg,  where  vg  is  defined  by  (3.29).  Then 
*  00 

f  e  L  (D;r) ,  and  the  problem 

Lv  =  f  in  D 


*  i  3v  , 

v  >r  3  TT  >r 

D  N 


has  unique  solution  in  V^D)  which  belongs  to  C  (D)  (by  Lemma  4. A). 
Clearly,  v*  =  wq  -  vq  is  the  solution  of  the  problem.  Hence  w^  e  C°(D). 


Proposition  4.7.  Assume  is  a  solution  of  (4.4).  Let 

=  { (r ,z)  e  D|r  >  rn,  w  <  g  (r,H)}  u  ft 
q  0  q  q  1 


Q  •  E»  D. 


ft  =  t(r,z)  e  D|r  >  rn,  w  >  g  (r,H) } 
q  0  q  q 


Then,  in  the  sense  of  distributions, 


-1  in  ft 


Lw  = 

q 


in  ft 


(4.23) 


*  00 

Proof:  By  (4.22)  both  ft  and  ft  are  open  sets.  Given  ♦  e  C.(ft 
-  q  q  0  q 

define  $  =  0  in  D\ft  .  Clearly  E.  C  ft  ,  where  E,  is  the  support  set 

q  V  q  y 

t>.  Let 


of 
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«.t ... -  ?****pKTm 


i 


m  *  min  (g  (r,H)  -  w  ) 

V(r>ro}  q  q 


Then  m  >  0,  and  there  exists  X  >0  such  that  for  each  real  X  with 


| X|  <  X  we  have 


|Xi|>|  <  m 


w  +  Xty  <  g  (r,H)  in  ft  ^  (r  >  r„) 
q  q  q  o 


On  choosing  in  (4.4)  v  *  wq  +  XiJ>,  we  obtain 


0  <  /q  rVw  *^(X40drdz  -  X  r<|»  drdz 


X  r^w  *Vi|)  drdz  *  X  r<|>  drdz 


As  the  sign  of  X  is  arbitrary,  we  have 


/n  r?w  drdz  *  /fi  r^  drdz  -V  if*-C  C~(ft  ) 


i.e.  Lw„  =  -1  in  ft  .  Similarly,  we  obtain  Lw„  =0  in  ft  . 

q  q  ■*  q  q 


Q.E.O. 


Lemma  4.B.  Let  f2  =  Lv2,  where  v2  is  defined  by  (3.31).  If  v  e  R(L)  , 


6  =  vf^  r  drdz  +0 

Proof :  Assume  @  =  0.  Then  f^  X  v,  and  f  i  R(L)^  (by  Theorem 

* 

f?  <3  R(L).  It  means  that  there  exists  v  8  tl(r>)  such  that 

(  Lv  *  f 


(4.24) 


4.B).  So 


00  *  J  _ 

By  (3.3R),  f  e  L  <D;r).  Hence  v  8  C  (D)  (by  Theorem  3. A).  Let 


/2  -  v  .  Then  v  3  U(D)  O  C  (D),  and 


i 

f 

Me 
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Lv  -  0 


in  D 


o  on  r  u  r2 

tn  —  on  r  u  r  u  r 

r„  3  4  5 


— 1 
3n  T 


v. 

Hence  v  has  minimum  in 
principle);  nor  on  r 


“  0  . 

N 

D,  which  lies  on  3D;  but  not  on 
(by  remark  3.2).  It  must  be  on 


r 

6 

and 


(by  Hopf 


min  v  =  v(rfl,h)  *  in  —  . 

D  1 


But  -y-  ■  0  on  I*  ,  then  t—  *  0  at  the  point  (r„,h)  (as  v  e  C  (D) . 
on  o  on  o 

This  contradicts  the  Hopf  principle. 


Q.E.D. 


Theorem  4.9.  Assume  vq  is  defined  by  (3.29),  fq  *  Lvq.  wq  is  the 

solution  of  (4.4),  F  «  Lw .  Let  v  e  R(L) 

'4  4 

G(q)  “  ID  (F  “  i  ) vrdrdz  (4.25) 

then  the  following  two  assertions  are  equivalent: 

( 1 )  q  is  a  root  of  the  equation 

G(q)  »  0  (4.26) 

(2)  w_  e  V2(D)  n  C^D) .  (4.27) 

q 

Proof :  If  G(q)  =  0,  then 

/  (F_  -  f_) vrdrdz  «*  0 

0  5  5 

i.e.  (F_  -  f__ )  1  v  in  V  (D).  By  Lemma  4. A  we  have 

q  q 

L(w_  -  v_)  e  R(L) 

q  q 

Hence 

w_  -  v_  e  U(D)  . 

q  q 
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It  follows  from  Theorem  3. A  that 


w__  -  v_  e  C^D)  .  (4.28) 

q  q 

We  obtain  (r.27).  Conversely,  if  (4.27)  is  valid,  then  F_  -  f_  e  R(L) • 

j_  q  q 

(4.26  follows  from  the  fact  that  v  e  R(L)  . 

Q.E.D. 

Lemma  4.10.  (4.26)  has  at  least  one  real  root. 


Proof :  At  first  we  prove  that  the  function 

F  (q)  *»  f  F  vrdrdz 
D  q 

is  continuous  in  -•  <  q  <  +“. 


Given  q*  e  R,  let  {q^  be  a  sequence  converging  to  q*,  wg  be  the 

* 

solution  of  (4.4)  corresponding  to  q* ,  v  be  defined  by  (3.29),  and  v 

x  qi  qi 


-  v  .  Then  we  have 
qi  qi 


30 


By  (3.29)  -  (3.32)  and  (4.11)  the  sequence  {f  -  f  }  is  bounded  in 


IT(Dir).  Therefore  it  is  possible  to  select  a  subsequence,  still  called 
(F  -  f  },  in  such  a  way  that 


It  follows  from  Lemma  4. A  and  Remark  4.4  that 


* 

v 


,  *  ,  *  1 
Iv  I  converges  strongly  to  v  in  V  (D) 


in  L2(Djr)  . 

(4.29) 

n  — 

in  C  (D) 

(4.30) 

in  (D) 

(4.31) 

where  v  satisfies 


Lv 


* 

v  I 


F 


3v 

D  N 


r_  ■  3n  'r  *  0 


Since  wq^  »  vq^  +  Vq  '  we  ^ave  (4.30)  and  (4.31)) 


(w  } 

converges  strongly  to  w 

in 

C°(D) 

(4.32) 

qi 

(w  } 

converges  strongly  to  w 

in 

V1(D) 

(4.33) 

qi 

* 

where  w  =■  vq,  +  v  . 

Accordingly,  w|p  *  g  , 
n  ^ 

(by 

(4.32)),  i.e. 

w  e  Kq,. 

We  also  have 

Lw  -  f  ,  +  F  .  (4.34) 

q' 

Now  fix  any  v  e  Kq,j  it  is  easily  seen  that  there  exists  a  sequence  {v^} 

such  that  Vj  e  K_  and  (v. }  converges  strongly  to  v  in  V^D);  so  that 
1  ”i  i 


from 


f  rVw  »7(v.  -  w  )drdz  -  (h  -  h)  J  (v  -  w  )| 
D  q,  i  q,  w  0  i  q.  1 


qA  I z=h 


rdr 


-  /D  r(vj  -  w'  )drdz  >  0 


-31- 


it  follows  (remark  that  (4.33)  implies  w'  ♦  w'  in  L* (Djr ) )  that 

gi 

r 

/  rVw»7(v-w)drdz  +  (h^-h)  /Q  (v-w)lz_h  r<Sr 

-  /  rfv'-w* )drdz  >0  . 

D 

Hence  w  is  the  solution  of  (4.4)  for  q  ■  q'»  i.e.  w  **  wqt,  an<*  (by 


(4.34)) 


Lw 


F 

q 


F  ,  -  f  ,  +  F 

q’  q' 


(4.35) 


By  (4.29)  and  (4.35)  we  have  that 


2 

{F  }  converges  weakly  to  F  ,  in  L  (Djr)  . 

n  a 


This  means  that  for  any  {q  }  with  lira  q  *  q*  there  exists  a  subsequence, 

still  called  tq^,  such  that 

lim  /  F  vrdrdz  -  /  F  ,  vrdrdz  . 


q*q 


_  D  Q 

+n 9  4 


D  q' 


+  it 

Accordingly,  F  (q)  is  continuous.  Clearly,  F  (q)  is  also  bounded  (by 


(4.11)). 


(4.26)  may  be  rewritten  as 

F  (q)  -  Bq  -  a  *  0  (4.26*) 

where  B  =  /Q  f2vrdrdz  +  0  (by  Lemma  4.8),  a  "  /D  f., vrdrdz.  Clearly,  the 
right  side  of  (4.26* )  changes  its  sign  when  q  changes  from  -•  to  +•• 
Hence  (4.26)  has  at  least  one  real  root. 


Q.  E.  D» 


1  —  2 

We  call  the  solution  of  (4.4)  regular  if  w  e  C  (D)  O  V  (D).  It  follows 

q 

immediately  from  Lemma  4.10  and  Theorem  4.9  that  the  following  theorem  is 
valid. 

Theorem  4.11.  There  exists  at  least  one  q  e  R  such  that  the  solution  w_ 

q 


of  (4.4)  is  regular. 
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(4.36) 


Proposition  4.12.  If  w_  is  a  regular  solution  of  (4.4),  then 

q 

q  <  q0 

where  qQ  is  defined  by  (4.10). 


Proofs  Let  w^  be  the  solution  of  (4.4),  and 

* 

F  (q,z)  *  w  (r„,z)  -  w  (r„,h)  -  (h  -h)(z-h)  for  z  <  h 
1  q  0  q  0  w 


F  (q)  *  lira 
*  z+h-0 


F1(q,z) 


Given  q  >  q^.  Assume  is  regular  solution.  Then  it  follows  by 

(4.15)  and  w  e  c'(D)  that 

q 

F2(q)  =*  0  .  (4.37) 

On  the  other  hand,  we  have 


Hence  w^  has  minimum  in  D,  which  lies  on  3d.  Clearly,  it  just  is 

w  (r  ,h)  <  0.  So  we  have 
q  0 

F^(q,z)  >  -(h^  -  h)(z  -  h) 

and 

F„(q)  <  -(h  -  h)  <  0  for  q  >  qA  .  (4.38) 

i  w  0 

This  contradicts  (4.37),  and  w^  cannot  be  a  regular  solution  for  q  >  q^. 

Q.E. D. 

To  complete  this  section  we  display  the  relation  between  (4.4),  (4.6)  and 
(4.8). 
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Theorem  4.13.  If  q  <  q^,  and  w^  is  the  solution  of  (4.4),  then 


(4.39) 


Moreover,  (4.4),  (4.6)  and  (4.8)  are  equivalent  for  q  <  q  . 

Proofs  By  Proposition  4.2  and  4.3  it  is  sufficient  to  prove  (4.39).  Let 
q  <  qQ  and  be  the  solution  of  (4.4).  It  follows  from  (4.11),  (4.15) 

(3.15)  and  (3.16)  that 


Lw  <  0,  w  | r 

q  qr„ 


0, 


3w 

_a, 

3n'r 


Ip  >  0 


N 


(4.40) 


It  is  easily  shown  by  the  maximum  principle  that  w^  >  0  in  D.  We  now  prove 

the  second  part  of  (4.39). 

* 

Let  ft^  be  defined  as  in  Proposition  4.7.  Noting  that 


and  w  e  C°(D) 

q 

and  Stampacchia 


w 

q 


w  in  ft. 

q  1 

max(w  ,g  (r,H)  in  D\ft. 
q  q  l 

S 

! 


we  may  prove  that  w’  e  v\d)  (see  for  instance  Kinderlehrer 

q 

[1980,  p.  50}).  Hence  w"  e  v'(D)  ^  c°(D).  We  have 

q 


Lw"  -  Lw  -  Lg  (r,H)  -  0  in  ft  (by  (4.2)  and  (4.23)) 

q  q  q  q 


w"|  .  “  0  (by  that  w"  e  C®(0)  and  w"|a_  *  0)  . 

ql3Q*  q  a  30 

q 

* 

It  follows  from  the  maximum  principle  that  w"  =  0  in  ft  .  Hence  w"  =  0  in 

q  q  q 

D  and  w_  *  w’ . 

q  q 

Q«E»  D* 
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It  follows 
Corollary  4.14. 


from  the  theorem  that 

Under  the  same  assumption  as  in  Theorem  4.13  we  have 


3w 


Ip  >  o 


(4.41) 


5.  The  Existence  of  the  solution  of  (PPW) . 

In  this  section  we  prove  that  a  regular  solution  of  (4.4)  corresponds  to 
a  solution  of  (PPW).  Following  the  framework  of  Baiocchi  et  al.  [1973],  we 
establish  several  lemmas  at  first. 

Throughout  this  section  let  wg  be  a  regular  solution  of  (4.4)  and  let 
ft  be  defined  as  in  Proposition  4.7. 

q 


Lemma  5.1. 


3w  _ 

^  0  in  D. 


(5.1) 


3w 


Proof;  Let  E  *  {(r,z)  e  D  0).  Then  E  is  an  open  set,  and  E  C  ft^. 

In  fact,  if  (r*,z*)  e  D  ft  ,  then  it  follows  from  (4.36)  and  (4.39)  that 

q 

*  *  *  * 

w  (r  ,z  )  =  g  (r  ,H)  >  w  (r  ,z)  for  0  <  z  <  H  . 

q  q  q 

_  *  * 

3w  ( r  ,  z  ) 


*  *  _ 


Hence 


0,  and  (r  ,z  )  e  E. 


If  E  +  0,  then  by  Proposition  4.7  we  have 


3w 

L^3l3^  =  0  in 


3w 


l.r. 


Therefore  t— ^  has  a  strictly  negative  minimum  in  E  (since  w  e  C  (D)) 
3z  3w  g 


which  lies  on  3E;  but  neither  on  3e  n  ft  where  *  0>  nor  on  (by 

3w  q 

(4.41));  nor  on  where  ■  h^  -  h  >  0;  nor  on  (by  remark  3.2). 


This  is  absurd,  and  E  -  0. 


Q.  E.  D. 
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Lemma  5.2.  If  q  >  0,  then 


3w 

0  <  <  SL  in  d  . 

3r  r 


(5.2) 


3w 

Proof.  Let  v  =  r  IT2*  ***  E  =  {(r,z)|v  >  q}.  Then  E  is  an  open  set,  and 
E  ft  (since  it  is  easy  to  show  that  v  <  q  in  D\ft  )•  Therefore, 

<1  g 


where  L1 


L  v  =  r(Lw  )  =0  in  E 

1  q  r 

2  3 

L  -  —  is  still  an  elliptic  operator. 


Simple  computation 


indicates : 


0 

v  * 

q 


3v 

3n  “  0  on 


on  ri  u  r? 
on  r3 


r  u  r  u 

2  4 


r5U 


By  maximum  principle  (see  for  instance  Gilbarg  and  Trundinger  [1977])  v  has 
maximum  strictly  bigger  than  q  in  E  which  lies  on  3e.  An  argument 
similar  to  that  in  the  proof  of  Lemma  5.1  indicates  that  E  *  0.  Similarly  we 
may  prove  that  E  =  {(r,z)|v  <  o)  =  0. 


Remark  5.1.  Similarly  we  have  that  if  q  <  0  then 

*  t—  <  0  in  D  . 
r  3r 


Remark  5.2.  If  q  ■  0,  then 


3wQ 

3r 


Hz 


2  2  (h  -z) 

z  H  w 


0.  Hence 

2 


for  h  <  z  <  h 


Q.  E . D . 


(5.3) 


It  requires  that  H  =  h  ..  We  have  assumed  that  h  <  H.  So  w„  is  not 

w  v/  0 

regular  solution. 
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y 


Remark  5.3.  By  using  (5.1),  (5.2)  and  (5.3)  we  may  easily  show  that  if 

(r/z)  e  D\n  then 
<3 

r  >  r„,  z  >  h  . 

0  w 

3w 

Lemma  5.3.  t—2  ■  0  on  T_  .  (5.4) 

“ •  oz  3 

Proofs  It  is  obvious  that 


3wq(ri«H) 

3z  ”  3z 


0 


(5.5) 


*  * 

Let  q  >  0.  Then  for  r  e  [r0»r^[  and  X  >  0  with  r  +  X  <  there 
exists  9  e  ]0,1[  such  that 


3w  (r  +X,h)  3w  (r  fH) 

- v— 1/x 


*  lim 
h++0 


.  w  (r  +X,H-h)-w  (r  +X,h) 

I  ,_2 _  _g_ _ 

X  1  -h 


*  * 
w  (r  ,H-h)-w  (r  ,H) 


-h 


*  -  lim 
h-*-+0 


1  T^w(r  +9X,H-h) 
h  1  3r 


-■  ■]  >  0  (by  Lemma  5.2)  . 

r  +9X 


Therefore  yj  is  nondecreasing  function  of  r  on  and  (5.4)  follows 

from  (5.5).  The  proof  for  q  <  0  is  similar. 

«  *  # 

For  point  p  =  (r  ,z  )  we  define  the  sets 

Q+«  =  {(r,z)  e  D|r  <  r  ,  z  >  z  } 

P 

5  *  =  t(r,z)  e  D|r  >  r  ,  z  <  z  }  . 

P 


Lemma  5.4.  If  q  >  0,  then 


Q+  C  d\8  for  p  e  D\i! 

p  q  q 


Q  C  8  for  p  6  D  n 
P  q  q 


(5.6) 


(5.7) 


*  —v  —  *  *  . 

Proofs  Let  p  e  D\8  and  a(r)  =  w(r,z  )  -  g(r,H).  We  have  r  >  rn, 

-  q  u 

z  >  h  (by  Remark  5.3),  o(r  )  =  0  (Theorem  4.1*),  o*(r)  <  0  (by  Lemma 

w 

* 

5.2).  Hence  o(r)  >  0  for  r  e  [rQ,r  ],  and  w(r,z)  *»  g(r,H)  for 
* 

r  e  (r  .r  1.  It  follows  from  Lemma  5.1  and  Theorem  4.13  that 
0 

*  *  i 

w(r,z)  =  g(r,H)  in  l(r,z)|r  *  r  <  r  ,  z  <z<H> 
i.e.  Q+.  C  d\8  . 

*  q 

P  ^ 

For  p'  *  (r*,z’)  e  D\Sg  there  exists  p  e  d\8^  such  that  r  >  r', 

z*  <  z'.  Clearly  Q+,  CQ+t  and  Q+*  c  D\8^.  Hence  Q+,  C  d\8^. 

P  P  P 

For  p  e  3d  (5.6)  is  trivial.  (5.7)  is  easily  seen  by  reducing  to 
absurdity  and  (5.6). 

Qe  Re  De 

*  *  * 

Remark  5.4.  For  q  <  0  and  p  *  (r  ,z  )  we  define 

R  *  -  l(r,z)  e  D|r  >  r  ,  z  >  z  > 

P 

_  *  * 

R  #  *  Ur,z)  e  D|r  <  r  ,  z  <  z  J  . 

P 

Then  we  obtain  by  similar  argument  that 

R+.  c  d\8  for  p  e  D\8 

p*  q  q 

R  .  r  8  for  p  e  D  n  9B  . 

p*  q  q 

From  Lemma  5.4  immediately  follows  a  property  of  8^. 

Corollary  5.5.  8^  is  a  connected  set. 

Lemma  5.6.  38  n  D  does  not  contain  any  vertical  or  horizontal  line  segment, 

-  q 

and  38  n  I*  *  0. 

q  3 
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Proof:  Assume  that  3fl  n  d  contains  a  vertical  line  seament  T  ■ 

-  q 

{(r,z)  e  D|r  ■  r',  z*  <  z  <  zH).  Denote 
=*  {(r,z)  e  D|r  >  r',  z'  <  z  <  z"},  N2  ■  {(r,z)  e  D|r  <  r',  z*  <  z  <  z"}.  Then 
K,  C  1)  and  N  C  d\0  (by  Lemma  5.4).  Hence  w  (r,z)  ■  g  (r,H)  and 

^  4  CJ  t  4 


3W  _ 

•  q/r  in  N2j  DWg 
the  Cauchy  problem 


-1  in  Nv 


Therefore 


is  the  solution  of 


r 


Lw  «  -1 

q 


in  N. 


wq'r*  “  9<r'*H> 


_-5L|  .St¬ 
ar  1  „*  r' 


By  the  uniqueness  of  the  solution  we  have 


Wq  =  ^2~  *n  +  4(r'2  “  +  9q(*,H)  in  ^  . 

2 

But  w  »  Hz  -  r—  on  r  .  This  contradiction  proves  that  an  n  D  does  not 
q  2  2  i 

contain  any  vertical  line  segment.  By  similar  argument  and  (5.4)  we  obtain 
that  30^  n  d  does  not  contain  any  horizontal  line  segment  and 

an  nr  -  0. 

q  3 

Q.E.D. 

Theorem  5.7.  If  q  <  0,  and 


n  -  n 

U  {(r,z)  e  D | r  >  r  ,  w  <  g  (r,H)} 

(S.8) 

q  i 

q  q  q 

*  (r)  » 

q 

8up{z|(r,z)  e  «q>  for  r  e  ]rQ,r1I 

(5.9) 

f  (r.)  3 

lim  <fi  (r),  <fi  <r„)  *  lim  <p  (r) 

(5.10) 

q  o 

r10.0  q  q  '  rt,-0  q 
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(5.11) 


3w 


u  ”  q— *  +  z  in  D,  u  ■  u  |n 
q  3z  q  q  « 


then  tu  ,  <fi  (r)}  is  the  solution  of  (PPW). 

q  q 

Proof :  First  we  note  that  ^(r)  is  a  well-defined,  strictly  increasing, 
continuous  function  for  r  6  Jr^r^l.  In  fact,  for  any  r  @  )r^,r^(  we  have 
wg  <  g^(r,H)  if  z  is  small  enough  since  wq  =  0  on  and  qq(r,H)  >0. 

So  (z|(r,z)  e  A  )  is  nonempty  and  (r)  is  well-defined.  It  immediately 

q  q 

follows  from  (5.1)  and  the  definition  of  ^g(r)  that 


fl  «  fl  u  {(r,z)  e  D|r  >  r.,  0  <  z  <  p  (r)}  . 

q  1  0  q 


(5.12) 


Lemma  5.6  shows  that  {(r,<0  (r))|r„  <  r  <  r.)  is  a  Lipschitz  graph  with 

q  0  i 

respect  to  the  axes  x  “  r-z,  y  *  r+z.  Hence  ipq(r)  is  a  strictly 

increasing,  continuous  function. 

By  virtue  of  (5.9),  (5.10)  and  (3.15)  it  is  readily  shown  that 

yp  (r„)  >  h  ,  V>(r4)  *  H.  Then  (2.6)  and  (2.7)  have  been  proved.  (2.8)  is 
Owi 

obvious . 

Now  we  check  (2.4).  Since  w_  =  in  D/S  we  have 

q  q  q 


3w 
_ c 

3z 


0,  u  =  z  on  T 

q  o 


the  rest  of  (2.4)  is  obvious  thanks  to  (3.15)  and  (3.16). 

2  — 

Finally  we  check  (2.9).  Given  any  ♦  e  C  (8  )  with  ^  -  0  in  a 

q 

neighborhood  of  U  (T^  n  9fl^)  U  u  we  have  (note  (4.23)  and  (5.12)) 


q 
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••Hit 


t*  [h (r  *f  !?  “  I;  ^r  a?  tlPldrdz 

q 

f  !!a  £i  „  \  H  „ 

“  ha  rJ^Tdrdr  +  r~drdidz 

q 

-  Ir^  r  ^  ||  dr  +  r  dz  {aince  -  0  on  U  r?) 

3  yf 

-q  fp  j-j-  dr  +  ^  dz  (since  w^  =  g^(r,H),  iyr^  *  q  in  D\0) 


-q['('<r0,v>(r0))  -  ♦(r1  ))  +  i|i(r0^(r0))  -  ,sp(T^  ))J  -  0 


The  proof  is  completed  by  virtue  of  the  denseness  of  (i|»)  described  above 


in  Kr 


Q.E.D. 


If  q  <  0,  then  by  using  similar  arguments  we  obtain  that  ^q(r)  is 
strictly  decreasing,  continuous  for  r  e  Jr  ,r  [,  and  that  >p  (r  )  ■  H,  This 


0'  1 

is  absurd.  Hence  we  obtain  (recall  (4.36)  and  remark  5.2): 
Proposition  5.8.  If  wg  is  a  regular  solution,  then 

0  <  q  <  qQ  . 

Proposition  5.9.  Let 


q'  V 


Then: 


Q  *  (q| w  is  a  regular  solution  of  (4.4)} 

q 


Q  c  10»q0i 


(5.13) 


(5.14) 


Q  is  a  closed  set 


Wq  is  nonincreasing  on  Q 


(5.15) 

(5.16) 

(5.17) 

Proof :  (5.15)  is  clear  by  virtue  of  (5.13).  (5.16)  follows  immediately  from 

Theorem  4.9  and  the  continuity  of  G(q)  (see  the  proof  of  Lemma  4.10).  Now 
we  prove  (5.17).  Let  q1#  q2  @  Q*,  q1  <  q2,  an<i 


E  *  {(r,z)  e  D|w 


<  O)  . 


Then  w  <  w  <  g  (r,H)  <  g  (r,H)  in  EH  (r  >  r  ).  Hence  ECO  ,  and 
q1  q2  q2  q1  0  q1 
Lw  <  0  in  E  (by  Proposition  4.4  and  4.7).  w  has  strictly  negative  minimum 

on  E  which  lies  on  3ej  but  not  on  3e  n  o  where  w  •  0»  nor  on  r 

where  w  >  Oj  nor  on  where  ^  ■  0»  nor  on  (by  remark  3.2).  This 


is  absurd.  Hence  E  *  0. 


Let  q  ■  igf(q),  q^  »  syp(q}.  By  (5.15)  and  (5.16)  we  have 
Q  ”  Q 


QeE« D  e 


%'<kt  e  0  '  >  °'  qM  <  q0 


From  (5.17)  follows  immediately  the  theorem 
Theorem  5.10.  For  any  q  e  Q*  we  have 


w  >w>w  inD. 

%  q  % 


(5.18) 


APPENDIX 


The  pi oof  of  Theorem  3. A 

Taken  6  with  0  <  6  <  rQ.  Let 

=  t(r,z)  e  D|0  <  r  <  6} 

<“2  »  {(r,z)  e  D|r  >  6/2} 

then  D  *  U  (o^.  Assume  P^,  is  the  corresponding  partition  of  unity 
then 

P1  +  P2  =  1  for  (r,z)  e  D  . 

We  may  choose  P_  such  that 


It  is  to  show  that 


L^v)  -  fP1  +  vLP1  +  27p1*7y 

g 

P1VI  -  ”a7(P1V)l  -  “  P1V'r=6  *  0  • 

1  >r  nu  3n  1  'r  nu  1 


It  follows  from  v  6  V^(D)  that  Vv  e  v^(w^).  Hence 

Vv  e  w’t^) 

*v 

where  is  three-dimensional  domain  obtained  by  rotating  around  z- 

axis.  By  the  embedding  theorem  we  have 


Vv  e  )  (p  <  6) 


So  Vv  e  Lp(m  >r)  and 


L(pfv)  e  lr(o>jir)  (p  <  6)  . 

Considering  P^v  as  the  solution  of  three-dimensional  problem  (2),  we  have 

p.v  e  W2{ u  )  (p  <  6)  . 

1  p  1 

By  using  the  embedding  theorem  again  and  returning  to  two-dimensional  dmain  we 


obtain 


P,v  e  C8(ui)  (0  <  |)  . 


In  o>2  *  polygon  A'BCEPQ'  the  operator  L  is  non-singular.  It  is 


easily  seen  that 


A(p2v)  =  g 


3(P2v) 

Vlrnv --K-Irn;  '  Vl  .  '  0 

D  2  N  2  2 


i\ 


where 


1  3  ^1 

g  -  P2f  +  vlp2  +  2Vp2«Vv  -  -  3^  (P2v)  e  L  (u>2)  (P1  <  6)  . 


Let  vo  *  9  *  w^ere  *  expresses  convolution  operation, 

/  r 


then  AvQ  -  g  and 


!  V 
.* 


*■  ,r  -Vw. •»*>£„;  ■  ••  Jfr.viy-**- a**  *$?-•, 


I 


NT! 


v0  e  w*  («2)  c  ck(u2)  (k  <  |)  . 


Let  v1  =  vQ  -  P2v,  then 


Av1  -  0 


V1*3w^pQ’  "  V0*3w^pQ‘ 


V1  vo 

3n~^pQ*  “  3n~lpQ*  * 


By  assumption  v  e  VZ(D)  we  have 


P2v  e  w2<<«2)  C  ck,(a)2>  (0  <  k*  <  ^)  . 


Hence 


v1  e  Ck’(a»2)  (0  <  k'  <  |)  . 

By  using  theorem  4.4  of  Volkov  (1965]  (Trudy  of  Mathematics  Institut  of 
Steklov,  77  (1965),  113-142)  we  have 


V1  e  c6(V  <  |>  • 


P2v  e  CB(«2)  <B  <  |)  . 


It  follows  from  (3)  and  (4)  that 


v  e  CB(D)  (B  <  |)  . 


QeEeDe 


A 
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